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Abstract. 

We compare V2/£2 and v^/e^ from single-shot and event-by-event (2+l)-dimensional hydrodynamic calculations and 
discuss the validity of using single-shot calculations as substitutes for event-by-event calculations. Further we present a proof- 
of-concept calculation demonstrating that v>2 and V3 together can be used to strongly reduce initial condition ambiguities. 

Keywords: Event-by-event, hydrodynamics, heavy-ion collisions, eccentricity, specific shear viscosity, elliptic flow, triangular flow 
PACS: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 

INTRODUCTION 

The Quark-Gluon-Plasma (QGP) created in heavy-ion collisions has been under intense study. In particular, it has 
been shown to exhibit almost perfect liquid collective behaviour. The extraction of one of its transport coefficients, the 
specific shear viscosity tj /s, has recently become one of the hottest topics. 

Shear viscosity reduces the conversion efficiency from initial geometry deformation to final flow anisotropies, and 
in [1] we reported that the eccentricity-scaled second and third order harmonic flow coefficients, V2/E2 and V3/E3, for 
unidentified charged hadrons are good choices to extract rj /s. For V2/62 this has already been done by several groups, 
but it was shown that this leaves a large uncertainty for rj/s due to ambiguities in the initial fireball deformation (see 
e.g. [2, 3]). We here show that a simultaneous analysis of V2/£2 and V3/E3 can resolve this ambiguity. The authors of 
[2, 3] use a "single-shot" hydrodynamic approach which evolves a smooth initial profile obtained by averaging over an 
ensemble of fluctuating bumpy initial conditions. We here address the question if it matters whether one instead follows 
nature's example and evolves each bumpy initial condition separately ("event-by-event hydrodynamics"), averaging 
over the fluctuating event ensemble only at the end. Due to the high numerical cost of the event-by-event approach 
on the one hand and the strong sensitivity of rj /s on V2/ £2 and V3 / £3 on the other hand this issue is of high practical 
relevance. 



SINGLE-SHOT VS. EVENT-BY-EVENT HYDRODYNAMIC CALCULATIONS 

We initialize the hydrodynamic simulations by generating fluctuating initial entropy density profiles from the Monte 
Carlo Glauber (MC-Glb.) and Monte Carlo KLN (MC-KLN) models for Au+Au collisions at yfs = 200 A GeV (see 
[1] for details). In Fig. 1 we compare the eccentricity-scaled elliptic and triangular flow coefficients for thermal pions, 
kaons and protons (i.e. without resonance decay contributions) from single-shot and event-by-event hydrodynamics, 
for ideal and viscous fluids. For the single-shot simulations, we average over the event ensemble in the participant 
plane, i.e. after centering and rotating each event by the participant plane angle ^23 ( e ) between the short axis of the 
second resp. third order harmonic component of its energy density profile e(x,y) and the impact parameter b [ 1 ] , before 
starting the evolution. 

For an ideal fluid (Fig. l(a,c)) event-by-event hydrodynamics produces significantly less elliptic and slightly less 
triangular flow than the equivalent single-shot evolution. The difference in V2/E2 is smallest for pions { ff(5%)) but 
increases with hadron mass to about 10% for protons. For v>3/£3 the differences are smaller, but again increase with 
hadron mass. In contrast, viscous hydrodynamics with 77 /s = 0.2 produces almost the same V2 / £2 and V3 / £3 for single- 
shot and event-by-event hydro, irrespective of hadron mass; if anything, the flows are now a little larger when the 
fluctuating events are evolved individually. Apparently, shear viscosity quickly damps the initial density fluctuations 
into something approaching the smooth ensemble-averaged profile before most of the flow develops. (The strange 
pattern seen at very large impact parameters in Figs. l(b,d) signals a jump of the flow angle xj/^ [1] by Tl/n (n = 2,3.) 
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FIGURE 1. Eccentricity-scaled elliptic {V2/E2, top) and triangular flow (v^/ei, bottom) as a function of impact parameter, for 
thermal pions, kaons, and protons from 200AGeV Au+Au collisions. Solid and dashed lines show results from single-shot and 
event-by-event hydrodynamics, respectively, for MC-KLN initial conditions. The left panels (a,c) assume ideal [1], the right panels 
(b,d) viscous fluid dynamic evolution with r\ fs = 0.2. 



Fig. 1 suggests that for a viscous fluid with sufficiently large viscosity rj/s single-shot hydrodynamics can substitute 
for event-by-event evolution for the calculation of both V2/e2 and V3 / £3 for unidentified and identified hadrons. For an 
ideal fluid this remains true for the triangular flow of unidentified charged hadrons (which are mostly pions) but not for 
that of identified heavy hadrons (e.g. protons), nor for the elliptic flow v^j&i- Until the QGP viscosity is known, it is 
therefore advisable to extract it from event-by-event hydrodynamic simulations. If it turns out large enough, additional 
systematic studies can be done using the more economic single-shot approach. 

By comparing the left and right panels in Fig. 1 we see that rj/s = 0.2 suppresses V3/e3 much more strongly (by 
~ 50%) than V2/£2 (which is suppressed only by ~ 25%). Taken together, V2/fi2 and V3/g3 thus over-constrain rj/s 
for a given model of initial state eccentricities £„. We will now show how this allows to distinguish experimentally 
between the MC-Glauber and MC-KLN models. 



REDUCING THE INITIAL CONDITION MODEL AMBIGUITY 

We reported in [1] that the MC-Glauber and MC-KLN models have similar £3 but the MC-KLN model has ~ 20% 
larger £2, and that for ideal fluid dynamics they give similar V2/£2 and V3/£3 ratios. Accordingly, for ideal fluids the 
two models generate similar V3 but different v%. We expect the ratios V2/fi2 and v^/ei, for these two models to remain 
similar even when adding viscosity (a corresponding study is ongoing). Therefore, for any fixed rj/s, they will generate 
similar V3 but different v%. Alternatively, when using different rj/s for the two models so that they produce the same 
V2, they will necessarily generate different V3. They will therefore not be able to simultaneously describe a given set 
of experimental V2 and V3 data with the same medium properties. 



To verify this statement quantitatively requires a proper event-by-event hydrodynmical calculation which is in 
progress. We here show instead results from a proof-of-concept calculation that supports our argument. We first 
generate a large set of deformed Gaussian initial conditions similar to the ones described in [4], but having both 
non-zero £2 and £3, as calculated from the MC-Glauber and MC-KLN models for the 20-30% centrality bin, with 
random relative orientation yf^ p — y/f p . We call these initial conditions "MC-Glauber-like" and "MC-KLN-like". 
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FIGURE 2. Differential V2(pr) (left) and V3 (pr) (right) from viscous hydrodynamics using MC-Glauber-like and MC-KLN-like 
initial conditions and different values for tj/s (see text for discussion). 

Fig. 2 shows differential V2.j(pt) curves resulting from the viscous hydrodynamic evolution of these initial condi- 
tions. The solid and dashed curves in the left panel show that, in order to obtain the same V2{pr) for MC-KLN-like and 
MC-Glauber like initial conditions, the fluid must be twice as viscous for the former than for the latter. The right panel 
shows that, with T] /s chosen to produce the same V2, MC-Glauber-like and MC-KLN-like initial conditions produce 
dramatically different V3, with the one from MC-KLN-like initialization being much smaller. Conversely, if rj/s is 
tuned to produce the same 1/3, MC-Glauber-like and MC-KLN-like initial conditions require the same value of rj/s 
(solid and dash-dotted lines in the right panel), which then leads to dramatically different V2 values for the different 
initial conditions (see corresponding lines in the left panel). These conclusions agree qualitatively with corresponding 
statements made in Refs. [5, 6]. 

SUMMARY 

We demonstrated that for sufficiently large viscosity {r\/s > ff(Q)2)) and limited precision requirements single-shot 
evolution of smooth averaged initial profiles can substitute for event-by-event evolution of fluctuating initial condi- 
tions, and that a simultaneous analysis of V2 and V3 overconstrains rj /s and thus has the power to discriminate between 
initial state models. A precise extraction of 77 /s without initial state ambiguity will, however, require event-by-event 
viscous hydrodynamical evolution of fluctuating initial states, coupled to a microscopic hadronic cascade for the 
freeze-out stage [3], to calculate ensemble averages of both vi and V3. 
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